EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2011-033 
2013/01/19 



CMS-SUS-10-004 



Search for new physics with same-sign isolated dilepton 
events with jets and missing transverse energy at the LHC 



The CMS CoUaboratioif] 



Abstract 



The results of searches for new physics in events with two same-sign isolated leptons, 
hadronic jets, and missing transverse energy in the final state are presented. The 
searches use an integrated luminosity of 35 pb^^ of pp collision data at a centre-of- 
mass energy of 7 TeV collected by the CMS experiment at the LHC. The observed 
numbers of events agree with the standard model predictions, and no evidence for 
new physics is found. To facilitate the interpretation of our data in a broader range 
of new physics scenarios, information on our event selection, detector response, and 
efficiencies is provided. 



Submitted to the Journal of High Energy Physics 



*See Appendix A for the list of collaboration members 



1 



1 Introduction 

Events with same-sign isolated lepton pairs from hadron collisions are very rare in the standard 
model (SM) but appear very naturally in many new physics scenarios. In particular, they have 
been proposed as signatures of supersymmetry (SUSY) IHHII, universal extra dimensions IH, 
pair production of T5 73 (a fermionic partner of the top quark) [Sj, heavy Majorana neutrinos HI, 
and same-sign top-pair resonances as predicted in theories with warped extra dimensions [Zl- 
In this paper we describe searches for new physics with same-sign isolated dileptons (ee, e}i, 
jiji, er, }iT, and rr), missing transverse energy {Ef^^^), and hadronic jets. Our choice of signal 
regions is driven by two simple observations. First, astrophysical evidence for dark matter [8j 
suggests that we concentrate on final states with E^'** . Second, observable new physics sig- 
nals with large cross sections are likely to be produced by strong interactions, and we thus 
expect significant hadronic activity in conjunction with the two same-sign leptons. Beyond 
these simple guiding principles, our searches are as independent of detailed features of new 
physics models as possible. The results are based on a data sample corresponding to an inte- 
grated luminosity of 35 pb^^ collected in pp collisions at a centre-of-mass energy of 7 TeV by 
the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) in 2010. 

This paper is organized as follows. The CMS detector is briefly described in Section 2] The re- 
construction of leptons, Ef^^^ , and jets at CMS is summarized in Section [s] Section |4 describes 
our search regions. We perform separate searches based on leptonic and hadronic triggers in 
order to cover a wider region in the parameter space of new physics. Electron and muon trig- 
gers allow for searches that require less hadronic energy in the event, while hadronic triggers 
allow inclusion of lower transverse momentum (pj) electrons and muons, as well as hadronic 
T decays in the final state. The dominant backgrounds for all three searches are estimated 
from data, as discussed in Section |5j Systematic uncertainties on the predicted number of sig- 
nal events and results of these searches are discussed in Sections |6] and [Zl We conclude with 
a discussion on how to use our results to constrain a wide variety of new physics models in 
Section H 



2 The CMS Detector 

A right-handed coordinate system is employed by the CMS experiment, with the origin at the 
nominal interaction point, the z-axis pointing to the centre of the LHC, and the y-axis pointing 
up (perpendicular to the LHC plane). The polar angle d is measured from the positive z-axis 
and the azimuthal angle (p is measured in the xy plane. The pseudorapidity is defined as rj = 
-ln[tan(f)]. 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, 13 m in length, providing an axial field of 3.8 T. Within the field volume are several par- 
ticle detection systems which each feature a cylindrical geometry, covering the full azimuthal 
range from Q < ^ < 2n. Silicon pixel and strip tracking detectors provide measurements of 
charged particle trajectories and extend to a pseudorapidity oi\r]\ = 2.5. A homogeneous crys- 
tal electromagnetic calorimeter (ECAL) and a sampling brass/ scintilator hadronic calorimeter 
(HCAL) surround the tracking volume and provide energy measurements of electrons, pho- 
tons, and hadronic jets up to \r]\ = 3.0. An iron-quartz fiber hadronic calorimeter, which is also 
part of the HCAL system, is located in the forward region defined by 3.0 < |//| < 5.0. Muons 
are measured in the pseudorapidity range oi \rj\ < 2.4, with detection planes made of three 
technologies: drift tubes, cathode strip chambers, and resistive plate chambers. These are in- 
strumented outside of the magnet coil within the steel return yoke. The CMS detector is nearly 
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hermetic, allowing for energy balance measurements in the plane transverse to the beam direc- 
tion. A two-tier trigger system is designed to select the most interesting pp collision events for 
use in physics analysis. A detailed description of the CMS detector can be found elsewhere [9l. 

3 Reconstruction of Leptons, Missing Energy, and Jets 

Muon candidates are required to be successfully reconstructed [lOJ using two algorithms, one 
in which tracks in the silicon detector are matched to consistent signals in the calorimeters and 
muon system, and another in which a simultaneous global fit is performed to hits in the silicon 
tracker and muon system. The track associated with the muon candidate is required to have a 
minimum number of hits in the silicon tracker, have a high-quality global fit including a min- 
imum number of hits in the muon detectors, and have calorimeter energy deposits consistent 
with originating from a minimum ionizing particle. 

Electron candidates are reconstructed fH} starting from a cluster of energy deposits in the 
ECAL, which is then matched to hits in the silicon tracker. A selection using electron identifica- 
tion variables based on shower shape and track-cluster matching is applied to the reconstructed 
candidates; the criteria are optimized in the context of the inclusive W — > ev measurement (121 
and are designed to maximally reject electron candidates from QCD multijet production while 
maintaining approximately 80% efficiency for electrons from the decay of W/Z bosons. Elec- 
tron candidates within AR = A(p^ + Atj^ < 0.1 of a muon are rejected to remove electron 
candidates due to muon bremsstrahlung and final-state radiation. Electron candidates origi- 
nating from photon conversions are suppressed by looking for a partner track and requiring 
no missing hits for the track fit in the inner layers of the tracking detectors. 

Hadronic r candidates (t;,) are identified IIT3l starting with a hadronic jet clustered from the 
particles reconstructed using the particle-flow global-event reconstruction algorithm |[T4L The 
highest-pr charged track within a cone of AR < 0.1 around the jet axis is required to have 
Pt > 5 GeV. A variable size cone of AR < 5 GeV/ pj is then defined around this track, and the 
boosted T-decay products are expected to be confined within this narrow cone. Only r candi- 
dates with one or three charged hadrons in this cone are selected. The discrimination between 
hadronic t decays and generic QCD jets is based on an ensemble of five neural networks, each 
of which has been trained to identify one of the five main hadronic T-decay modes using the 
kinematics of the reconstructed charged and neutral pions llTSl . 

All lepton candidates are required to have |?/| < 2.4, and be consistent with originating from 
the same interaction vertex. Charged leptons from the decay of W/Z bosons, as well as the 
new physics we are searching for, are expected to be isolated from other activity in the event. 
We calculate a relative measure of this isolation denoted as Rellso . This quantity is defined 
as the ratio of the scalar sum of transverse track momenta and transverse calorimeter energy 
deposits within a cone of AR < 0.3 around the lepton candidate direction at the origin, to the 
transverse momentum of the candidate. The contribution from the candidate itself is excluded. 

In order to suppress the background due to dileptons originating from the same jet, we require 
that selected dileptons have a minimum invariant mass of 5 GeV. This helps to keep dileptons 
uncorrelated with respect to their Rellso observables, which is a feature we exploit in the anal- 
ysis. We also remove events with a third lepton of opposite sign and same flavour as one of the 
two selected leptons if the invariant mass of the pair is between 76 and 106 GeV. This require- 
ment further reduces an already small background contribution from WZ and ZZ production. 

Jets and E^'^^are reconstructed based on the particle-flow technique desribed in l[T4llT6| . For 
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jet clustering, we use the anti-Zcr algorithm with the distance parameter R = 0.5 lllTll . Jets are 
required to pass standard quality requirements flS] to remove those consistent with calorimeter 
noise. Jet energies are corrected for residual nonuniformity and nonlinearity of the detector 
response derived using collision data I1T9| . We require jets to have transverse energy above 
30 GeV and to be within \tj\ < 2.5. We define the Hj observable as the scalar sum of the pj of 
all such jets with AR > 0.4 to the nearest lepton passing all our requirements. 



The searches discussed in this paper employ two different trigger strategies, electron and muon 
triggers in one case, and Hj triggers in the other. The leptonic triggers allow for lower Hj 
requirements, while the Hj triggers allow for lower lepton pj, as well as final states with 
hadronic t decays. The motivation for covering the widest possible phase space in this search 
can be illustrated by an example of a SUSY cascade, shown in Fig. [l| naturally giving rise to 
jets, and same-sign leptons: (gluinos/squarks) (charged gaugino) (lightest super- 

symmetric particle (LSP) neutralino). The mass difference between the gluino/ squarks and the 
charged gaugino, typically arbitrary, defines the amount of hadronic activity one may expect 
in the event. The mass difference between the gaugino and a neutralino influences the lepton 
Pt spectrum. Depending on the nature of the chargino and neutralino, their mass difference 
can be either arbitrary (e.g., wino and bino) or typically small (e.g., higgsinos). Moreover, there 
are a number of ways to generate a large production asymmetry between r and e/}i leptons, 
which motivates us to look specifically for events with a t. 



Figure 1: An example of a process involving the production and decays of SUSY particles, 
which gives rise to two same-sign prompt leptons, jets, and missing transverse energy. 

In the following we describe the search regions explored by each trigger strategy. As a new- 
physics reference point, we use LMO, a point in the constrained Minimal Supersymmetric Stan- 
dard Model (CMSSM) [20] defined with the model parameters mo = 200 GeV, mi/2 = 160 GeV, 
tanj6 = 10, ji > 0, and Aq = —400 GeV. LMO is one of the common CMSSM reference points 
used in CMS across many analyses. As the abbreviation suggests, LMO provides squarks and 
gluinos with relatively low masses, and thus has a large production cross section. It is beyond 
the exclusion reach of the searches performed by LEP and Tevatron, but it has recently been 
excluded by searches pT] - t23l with ATLAS and CMS concurrently with this one. Nonetheless, 
we continue to use LMO as it provides a common model for which to compare our sensitivity 
with that of other analyses. Aside from the LMO point, several SM simulation samples are used 
to both validate and complement various background estimation methods that are based on 
the data itself. These samples rely on either PYTHIA 6.4 [24] or MADGRAPH [25] for event 
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generation and GEANT4 l|26| for simulation of the CMS detector. Samples used include tt, 
single-t, 7+jets, W+jets, Z+jets, WW, WZ, ZZ, and QCD multijet production. Next-to-leading- 
order (NLO) cross sections are used for all samples except for QCD multijet production. 

4.1 Searches using Lepton Triggers 

We start with a baseline selection inspired by our published tt — ?► + X (£ = e or cross 
section measurement ||27|. 

Events are collected using single and dilepton triggers. The detailed implementation of these 
triggers evolved throughout the 2010 data-collecting period as the LHC instantaneous lumi- 
nosity was increasing. Trigger efficiencies are measured from a pure lepton sample collected 
using Z — > decays from data. The luminosity-averaged efficiency to trigger on events 

with two leptons with \tj\ < 2.4 and pj > 10 GeV, one of which also has pj > 20 GeV, is 
very high. For example, the trigger efficiency for an LMO event passing the baseline selection 
described below is estimated to be (99 ± 1)%. 

One of the electrons and muons must have pj > 20 GeV and the second one must have pj > 
10 GeV. Both leptons must be isolated. The isolation requirement is based on the Rellso vari- 
able introduced earlier. We require Rellso < 0.1 for leptons of pT > 20 GeV, and the isolation 
sum (i.e., the numerator of the Rellso expression) to be less than 2 GeV for pj < 20 GeV. 

We require the presence of at least two reconstructed jets, implying Hj > 60 GeV. Finally, we 
require the missing transverse energy £^^^^> 30 GeV (ee and or E^^^^> 20 GeV {e}i). This 
defines our baseline selection. 

Following the guiding principles discussed in the introduction, we define two search regions. 
The first has high E^'"" {E^'^^y 80 GeV); the second has high Ht {Hj > 200 GeV). These 
E^** and Hj values were chosen to obtain an SM background expectation in simulation of 1 / 3 
of an event in either of the two overlapping search regions. 
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Figure 2: Hj versus E^^^^ scatter plots for baseline region. (Left) Overlay of the three observed 
events with the expected signal distribution for LMO. The three observed events all scatter in 
the lower left corner of the plot. (Right) Scatter plot of the background in data when only one 
of the two leptons is required to be isolated. 
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Figure|2]shows the Hj versus £^®^ scatter plot for the baseUne selection, indicating the £^^^ and 
Ht requirements for the two search regions via horizontal and vertical lines, respectively. Fig- 
ure |2] (left) shows three events (red dots) in the baseline region, one of which barely satisfies 
Hj > 200 GeV, but fails the E™*® > 80 GeV requirement. In contrast, most of the signal from 
typical supersymmetry models tends to pass both of these requirements, as is visible in the 
LMO expected signal distribution overlaid in Fig. |2] (left). Backgrounds to this analysis are 
dominated by events with jets mimicking leptons, as discussed in Section |5j Requiring only 
one of the two leptons to be isolated thus allows us to increase the background statistics in 
order to display the expected distribution of SM background events in the (£^'^®, Ht) plane, as 
shown in Fig. (right). Backgrounds clearly cluster at low E^'^^and low Hj, with slightly more 
than half of the events failing both the E™^^ and Hj selections. Moreover, comparing the left 
and right plots in Fig. |2] indicates that the lepton isolation requirement on both leptons versus 
only one lepton reduces the backgrounds by roughly a factor of ten. 

4.2 Searches using Hadronic Triggers 

Hadronic triggers allow us to explore the phase space with low-p^ electrons and muons, as 
well as final states with hadronic T decays. We allow muons (electrons) with pj as low as 5 (10) 
GeV, and restrict ourselves to T;, with visible transverse momentum > 15 GeV, where T;, refers 
to hadronic t candidates only. All leptons must be isolated with Rellso < 0.15. 

For the ee, e}i, and final states, we require at least two jets, Ht > 300 GeV, and E^^^^ > 30 
GeV. As backgrounds from QCD multijet production are significant for t;,, we increase the 
£^^*^ and Hj requirements to £^'^^ > 50 GeV and Hj > 350 GeV in the ez;,, jiTh, and TjjTh final 
states. 

Figure |3] shows the efficiency turn-on curves for the Hj triggers used during three different 
data taking periods. The trigger thresholds were changing in order to cope with the increasing 
instantaneous luminosities over the 2010 running period. Roughly half of the integrated lumi- 
nosity in 2010 was taken with the highest threshold trigger. This measurement indicates that at 
Hj = 300 GeV the efficiency reaches (94 ±5)% . These trigger turn-on curves are measured in 
data with events selected by muon triggers. 
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Figure 3: Hj Trigger efficiency as a function of the reconstructed Hj for three data-collecting 
periods: 7 pb^^ with Hjl, 10 pb^^ with Hj2, and 18 pb^^ with Hj3. 
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5 Background Estimation 

Standard model sources of same-sign dilepton events with both leptons coming from a W or 
Z decay are very small in our data sample. Simulation-based predictions of the combined 
yields for qq — > WZ and ZZ, double "W-strahlung" qq — > q'q'W^W^, double parton scattering 
2 X (qq — ?■ W^), ttW, and WWW comprise no more than a few percent of the total background 
in any of the final states considered. As these processes have never been measured in proton- 
proton collisions, and their background contributions are very small, we evaluate them using 
simulation, assigning a 50% systematic uncertainty. The background contribution from pp — ?► 
W7, where the W decays leptonically and the photon converts in the detector material giving 
rise to an isolated electron, is also estimated from simulation and found to be negligible. All 
other backgrounds are evaluated from data, as discussed below. 

Backgrounds in all of our searches are dominated by one or two jets mimicking the lepton sig- 
nature. Such lepton candidates can be genuine leptons from heavy-flavour decays, electrons 
from unidentified photon conversions, muons from meson decays in flight, hadrons recon- 
structed as leptons, or jet fluctuations leading to hadronic t signatures. We will refer to all of 
these as "fake leptons". Leptons from W, Z, gauginos, etc., i.e., the signal we are searching for, 
will be referred to as "prompt leptons". 

The dominant background contribution is from events with one lepton, jets, and E^'^® — mostly 
tt with one lepton from the W decay, and a second lepton from the decay of a heavy-flavour 
particle. These events contain one prompt and one fake lepton, and are estimated via two 



different techniques described in Sections 5.1 and 5.2 While both techniques implement an 



extrapolation in lepton isolation, they differ in the assumptions made. Both techniques lead 
to consistent predictions as described in Section [5.4} providing additional confidence in the 
results. Backgrounds with two fake leptons are generally smaller, except in the final state with 
two hadronic r leptons, where the dominant background source is QCD multijet production. 
Contributions due to fake are estimated using an extrapolation from "loose" to "tight" t;. 



identification, as described in Section 5.3 



For the ee and eji final states, electron charge misreconstruction due to hard bremsstrahlung 
poses another potentially important background, as there are significant opposite-sign ee and 
e}i contributions, especially from tt, where both W's from the top quarks decay leptonically. 
This is discussed in Section |5^ 



5.1 Searches using Lepton Triggers 

Contributions from fake leptons are estimated using the so-called "tight-loose" (TL) method f27[ 
|28]. In this method the probability ctl for a lepton passing loose selections to also pass the tight 
analysis selections is measured in QCD multijet events as a function of lepton pj and rj. The 
key assumption of the method is that eji is approximately universal, i.e., it is the same for all 
jets in all event samples. Tests of the validity of this assumption are described below. 

The main difference between the tight and loose lepton selections is that the requirement on 
the Rellso variable defined in Section |3] is relaxed from Rellso < 0.1 to Rellso < 0.4. Other 
requirements that are relaxed are those on the distance of closest approach between the lepton 
track and the beamline (impact parameter) and, in the case of muons, the selection on the of 
the muon track fit. 

The quantity eji is measured in a sample of lepton-trigger events with at least one jet satis- 
fying Pt > 40 GeV and well separated (AR > 1) from the lepton candidate. We refer to this 
jet as the "away-jet". We reduce the impact of electroweak background (W, Z, tt) by excluding 
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events with 7. ^ U candidates, events with > 20 GeV, and events where the transverse 
mass Mr of the lepton and the E^^^ is greater than 25 GeV. Studies based on simulation indi- 
cate that this procedure results in an unbiased estimate of eji up to lepton pj ^ 40 GeV. At 
higher transverse lepton momenta the remaining electroweak contributions in the sample have 
a significant effect. Thus, eji is measured only up to pj = 35 GeV. It is taken to be constant at 
higher transverse momenta, as suggested by simulation studies. 
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Figure 4: Electron (left) and muon (right) TL probability ejL computed from QCD multijet 
events with different requirements on the minimum pj of the away-jet. The probabilities 
shown are projections of the two-dimensional function ejiif], Pt) onto the pj axis. 



The level of universality of Ctl is tested with different jet samples. Two types of tests are 
relevant and both involve the parent jet from which the lepton originates. The first test explores 
the sensitivity to the jet's pj, and the second test explores the sensitivity to the jet's heavy- 
flavour content. 

Sensitivity to jet pr stems from the fact that the probability for a lepton of a given pj to pass 
the Rellso selection depends on the pj of the parton from which the lepton originates. To be 
explicit, a 10 GeV lepton originating from a 60 GeV b quark is less likely to pass our Rellso re- 
quirement than the same lepton originating from a 20 GeV b quark. The heavy flavour sensi- 
tivity can be traced to semileptonic decays, which are a source of leptons in bottom and charm 
jets, but not in light-quark and gluon jets. 

To test the jet pr dependence, we select loose leptons in events with the away-jet above a vary- 
ing jet Pt threshold. Since these events are mostly QCD dijets, the pj of the away-jet is a good 
measure of the pr of the jet from which the lepton originates. We weight each event by eji mea- 
sured as described above, i.e., requiring the away-jet to have pr > 40 GeV. We then sum the 
weights and compare the sum to the number of observed leptons passing tight requirements. 
Varying the away-jet minimum pj requirement from 20 to 60 GeV, we find the observed yield 
to differ from the predicted yield by +54% (+49%) and —4% (—3%) for muons (electrons) in 
this test. The percentages here, as well as throughout this section refer to (observed — pre- 
dicted) /predicted. This non-negligible jet pr dependence can also be seen in Fig. |4| where we 
show €tl calculated using different away-jet thresholds. To test the heavy-flavour dependence, 
we repeat the exercise requiring that the away-jet be above pj > 40 GeV and be b-tagged, i.e., 
a jet in which we find a secondary vertex well separated from the interaction point consistent 
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with a b-hadron decay. By applying the b tag on the away-jet, the sample of jets from which 
the lepton originates is enriched in heavy flavours. Introducing this b tag we find the observed 
yield to differ from that predicted by —3% (—15%) for muons (electrons). We have thus shown 
that applying an eji obtained without a b-tagging requirement to a sample with such a re- 
quirement leads to a modest difference between observed and predicted. This validates our 
assumption that the TL method is flavour universal. 

To predict the background from prompt lepton + jets events in a signal region, the TL proba- 
bility is applied to a sample of dilepton events satisfying all the signal selection requirements, 
but where one of the leptons fails the tight selections and passes the loose ones. Each event is 
weighted by the factor ctl/ (1 — £tl), where eji is the tight-to-loose probability for the loose 
lepton in the event. The background contribution from this source is then estimated by sum- 
ming the weights of all such events (Si). 

The sum Si also includes the contribution from backgrounds with two fake leptons. However, 
these are double counted because in the case of two fake leptons passing loose requirements 
there are two combinations with one lepton passing the tight selections. The background con- 
tribution with two fake leptons is estimated separately by selecting events where both lep- 
tons pass the loose requirements but fail the tight requirements. Each event in this sample is 
weighted by the product of the two factors of cji / {1 — eji) corresponding to the two leptons in 
the event, and the sum ^2 of weights is used to estimate the background with two fake leptons. 

The total background from events with one or two fake leptons is then obtained as Si — Si- In 
kinematic regions of interest for this search, ^2 is typically more than one order of magnitude 
smaller than Si, indicating that the main background contribution is from one prompt lepton 
and one fake lepton. 

The method has been tested on simulated tt and W + jets events. In these tests, we use €tl 
measured from QCD simulation events to predict the number of same-sign dilepton events 
in these samples. In the tt simulation sample, we find that the observed yield differs from 
the prediction for the baseline selection by -41% (fake muons) and -47% (fake electrons). Ob- 
served and predicted yields in this test are consistent with each other at the 5% confidence 
level (CL) because the simulation statistics are modest. The same level of agreement is also 
found for the two search regions. In the W + jets case, the ratio of predicted events to ob- 
served is 0.8 ± 0.4 for fake electrons; the statistics for fake muons are not sufficient to draw 
any definitive conclusions. Based on these studies, as well as the dependence on away-jet pj 
and heavy-flavour composition discussed above, we assign a ±50% systematic uncertainty on 
the ratio (observed — predicted) /predicted and, hence, on the estimation of backgrounds due 
to fake leptons. In addition to this systematic uncertainty, the method has significant statisti- 
cal uncertainties based on the number of events in the samples to which the TL probability is 
applied. We find 6 (4) events in these samples for the £^'®* > 80 GeV (Hj > 200 GeV) search re- 
gions. The resulting background estimates in the two regions are 1.1 ± 0.6 and 0.9 ± 0.6 events, 
respectively, including only statistical uncertainties. 

As an additional cross-check, we determine the background estimate and observed yields in 
the baseline region. We estimate 3.2 ± 0.9 ±1.6 events from background due to fake leptons 
alone, and 3.4 ± 0.9 ±1.6 after all backgrounds are taken into account. The uncertainties here 
are statistical and systematic, respectively. The composition of the total background is esti- 
mated to be 86% (7%) events with one (two) fake leptons, 3% due to charge misidentification, 
and 4% irreducible background for which both leptons are isolated leptons from leptonic W 



or Z decay. As mentioned in Section 4.1 there are 3 events observed in the baseline region. 



in good agreement with the background estimate. Applying E^^^ > 80 GeV or Ht > 200 GeV 
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increases the fraction of events with one fake lepton, but statistical uncertainties on the individ- 
ual components of the background estimate are too large to meaningfully quantify the change 
in the relative contributions of the components. 

5.2 Search using Hadronic Triggers, Electrons, and Muons 



As described in Section 4.2 hadronic triggers allow us to explore the phase space with low-p^ 
leptons. However, lowering the lepton pj is expected to increase the relative contribution of 
events with two fake leptons. As shown below, this background now constitutes roughly 30% 
of the total background, as compared to only a few percent for the higher-pj thresholds in the 
search regions for the leptonic trigger analysis. This motivates the development of a method 
exclusively dedicated to predicting and understanding the QCD background with two fake 
leptons. 

At the same time, increasing the Hj requirement to 300 GeV, as driven by the hadronic trigger 
thresholds, reduces the expected W+jets background to only a few percent of the total back- 
ground. The background with one fake lepton is now reduced to tt and single-t processes, 
where the fake lepton is due mostly to semileptonic b decays. Therefore, we tailor the method 
for estimating the background with one fake lepton to these expectations. The method is sim- 
ilar to the TL technique, but has a number of important differences that are discussed further 
below. 

For this analysis, the estimation of the background with fake leptons starts with an evalua- 
tion of background events with two fake leptons and then proceeds with an estimation of the 
contribution of events with only one such lepton. 

First, we define a preselection control sample of events from the Hj-triggered data stream with 
at least two same-sign dileptons and with all event selection requirements applied, except for 
those related to E^'^^and isolation. We find 223 f/^, 6 ee, and 78 e^ events of this type. The large 
asymmetry between muons and electrons is mostly due to differences in the corresponding pj 
thresholds of 5 and 10 GeV, respectively. In addition, identification of electrons within jets is 
less efficient than that for muons. 

The preselection control sample is dominated by QCD multijet production. Studies based on 
simulation suggest that we should attribute about 10% of the preselection yields to tt contami- 
nation, while attributing a much smaller fraction to W+jets. 

The contribution from events with two fake leptons to the signal region is estimated by assum- 
ing that the three requirements, Rellso < 0.15 for each lepton and £^'*® > 30 GeV, are mutually 
independent and, hence, the total background-suppression efficiency can be written in the fac- 
torized form etot = e^jiso " ££2130 " ^met- This assumption has been verified both in simulation 
and directly in data. With simulation it is straightforward to prove the principle in the nominal 
preselection region because we can safely measure the efficiencies in a dedicated QCD sample 
where we know all leptons can be considered background (i.e., no contamination from prompt 
leptons exists). In data the contribution from prompt leptons is non-negligible and therefore 
some extra selection requirements are necessary to isolate a QCD enriched control sample. 

We validate the factorized expression for etot in two steps in data. First, we demonstrate that the 
selection requirement on Rellso is independent for each lepton. We begin by relaxing the Hj 
selection to 200 GeV and add events collected with leptonic triggers to gain more statistics. We 
then require £^'** < 20 GeV to suppress events with leptonic W decays. Figure [s] (left) shows 
that the single-muon efficiency can be squared to obtain the double-muon efficiency, thus vali- 
dating the assumption that the Rellso observable is uncorrelated between the two fake leptons 
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and the efficiencies can be factorized. In the second step, we demonstrate that the £™^^*and 
Rellso selection requirements are mutually independent. To accomplish this in data, we main- 
tain Ht above 300 GeV, but we include single-lepton events to increase statistics. To suppress 
the contributions from events with leptonic W decays, we modify the selection requirement on 
the lepton impact parameter from the nominal do < 0.2 mm to do > 0.1 mm. Figure |5] (right) 
shows that the Rellso selection efficiency for muons and electrons remains constant as a func- 
tion of the £^^^ selection requirement. The dashed lines represent the zeroth-order polynomial 
fits to the efficiency measurements made in the do control region for muons and electrons, 
respectively. For completeness, we also show the obvious bias arising when the impact param- 
eter requirement is inverted to do < 0.1 mm to enrich the sample with leptonic W decays. It 
is important to note that no attempt is made to apply the Rellso selection efficiency measured 
in the control region defined by do > 0.1 mm to the above formula for etot- This control region 
is only used to demonstrate the stability of the Rellso selection efficiency with respect to the 
£^^^^ requirement. The actual values of e^^iso and e^^iso ^re measured in the nominal preselec- 
tion region {do < 0.2 mm), where we assume this stability, and hence factorization, remains 
valid for events with two fake leptons. 

Having validated the selection factorization hypothesis, we proceed to measure the isolation 
and E^^^®^ selection efficiencies, one at a time, in the preselection control sample, where we ob- 
tain e,,iso = 0.036 ± 0.015, eeiso = 0.11 ± 0.08, eMET = 0.27 ± 0.03. Uncertainties quoted are 
statistical only. As before, we suppress leptonic W decays to reduce possible biases. We accom- 
plish this by requiring either E^^^^< 20 GeV or Rellso > 0.2 when measuring e^aso/ ^eiso/ or 
^MET- The appropriate product of these efficiencies is then applied to the event counts observed 
in the preselection control sample, leading to the background estimate of 0.18 ± 0.12 ± 0.12 
events. Uncertainties here are statistical and systematic, respectively. 
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Figure 5: (Left) The lepton isolation efficiency for one (solid squares) and two (open squares) 
leptons as a function of the relative isolation parameter cut. Also shown is the predicted 
double-lepton efficiency if the two lepton efficiencies are assumed to be independent of each 
other. Only the dimuon sample is shown here. (Right) The lepton isolation efficiency as a func- 
tion of the £™®®cut for electrons and muons with different requirements on the lepton impact 
parameter. Details are given in the text. 



The systematic uncertainties quoted above have two dominant sources. One is due to limited 
statistics in simulation and data when validating that the three requirements are indeed inde- 
pendent. We take the statistical precision (25%) of this cross-check as our systematic uncertainty 
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on the method. The other dominant source of systematic uncertainty can be attributed to the 
inabihty of the inverted selection requirements on E^'^^and Rellso to fully suppress contribu- 
tions from leptonic W decays (e.g., tt, W+jets) while measuring e^iso, ^eiso/ arid eMET- Studies 
based on simulation suggest that the bias (overestimate) can be as large as 60%, mostly via a 
bias in measuring iso- Conservatively, we do not correct for the possible bias, but take it as 
a systematic uncertainty. We thus arrive at a 65% systematic uncertainty on the estimate of 
backgrounds due to events with two fake leptons by adding these two effects in quadrature. In 
addition, we verified in simulation that the techniques used to suppress leptonic W decays (i.e., 
inversion of the E^'^^and Rellso requirements) do not alter or bias the selection efficiencies for 
fake leptons. 

It is worth mentioning that this method of evaluating background with two fake leptons does 
not require any reweighting of measured efficiencies. The average efficiencies are obtained 
from a QCD-dominated subset of the preselection sample, and then applied directly to the 
preselection sample as a whole to derive the prediction for the number of events with two fake 
leptons in the signal region. 

Next, we proceed with estimating the contribution of backgrounds with a single misidentified 
lepton. We start from a tight-loose control sample, to be further referred to as a sideband, and 
use the isolation selection efficiency for b jets, referred to as e*^*^), to predict event counts in the 
signal region. The sideband control sample contains events passing all signal selection criteria, 
except one of the two leptons is now required to have Rellso > 0.15. To begin, we count the 
number of events in this sample: 11 2 (ee), 6 (ef/), and 5 i}ie), the last lepton indicating 
which one in the pair is non-isolated. Then, we estimate the contribution of the background 
with two fake leptons to the sideband sample using the efficiencies quoted above from the fac- 
torization procedure. For example, for the dimuon channel, the contribution to the sideband 
from events with two fake leptons is N^;^ preselected ■ 2ef(iso(l — ^fdso) • ^met- The resulting yield 
estimates for events with two fake leptons are 4.2 (f/^w), 0.32 (ee), 2.3 (e/^), and 0.68 (fie). After 
subtracting this contribution, the remaining yields in the sideband are consistent with simula- 
tion predictions assuming that only tt (76%), single-t (7%), and W+jets (15%) contribute. This 
remaining sideband yield after subtraction is then scaled by an appropriate factor determined 
using the BTag-and-probe method |28|, as described below. 

The BTag-and-probe method relies on the basic premise that events with one fake lepton can 
be attributed to tt production, with one prompt lepton from leptonic W decay and the second 

fake lepton from semi-leptonic b decay. The efficiencies ej^^jg^ and e^'^g^ are thus defined as 
the probabilities of a muon or electron from semi-leptonic b decay to pass the Rellso < 0.15 
selection. These efficiencies can be measured in data using appropriately selected events from 

bb production. To determine e^^^^^ and ef^^g, we select a bb enriched control sample by requiring 
one b-tagged away-jet and one lepton candidate. In addition, we require Ht > 100 GeV to 
arrive at a b-quark pj spectrum similar to that expected for the tt background in the search. To 
reduce the bias from leptonic W or Z decays, we furthermore require E™'*® < 15 GeV and Mj < 
15 GeV, and veto events with two leptons forming a mass within 7 GeV of the mass of the Z 
boson. Approximately 80% of the leptons in this sample are from semileptonic heavy-flavour 
decay. 

We find that the resulting bb control sample differs sufficiently from the expected tt background 
in both lepton kinematics and jet multiplicity (Njets) to warrant corrections. We therefore mea- 
sure the Rellso distribution in the bb control sample in data in bins of lepton pj and Njets/ arid 
reweight these distributions using event probabilities co{pt, Njets) derived from a tt simulation 
sample. The resulting reweighted Rellso distributions for these three samples are overlaid in 
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Fig. |6] for muons (left) and electrons (right). The plots show distributions for tt simulation (red 
crosses) after all selections except Rellso on one of the two leptons, reweighted bb simulation 
(grey shade), and reweighted bb-enriched data (black dots). The agreement between the two 
simulation-based distributions validates the method. Agreement between data and simulation 
is observed but is not required for this method to be valid. The contents of the first bin of the 

two data plots are the above mentioned isolation selection efficiencies e^^i^^ and ^e'iso- 
e^^l^ = 0.029^°;°°^ and e^^^^ = 0.036+° with uncertainties due to statistics only. 
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Figure 6: Isolation variable distributions obtained with the BTag-and-probe method for muons 
(left) and electrons (right). Efficiencies for the Rellso < 0.15 (first bin in the distributions 
shown) are explicitly quoted. 

We probed four different potential sources of systematic uncertainties in the BTag-and-probe 
method, and added their contributions in quadrature to arrive at a total systematic uncertainty 

on the e^^^ efficiency of 54 (29)% for electrons (muons). The largest contribution to this uncer- 
tainty is due to the statistical precision with which the method in simulation is verified. The 
subdominant contributions include the change in the background prediction when the Hj re- 
quirement for selecting the bb control sample in data is varied from 100 to 150 GeV, or when 
the W-i-jets contribution is assumed twice the size predicted by the simulation with a value of 

e-^^ that differs from tt by a factor two, and the effect of reweighting the simulated events to 

lower the fraction of leptons not from b decay in the sample used to measure . We then mul- 
tiply the sideband yield in each of the four channels fifi, ee, e}i, and }ie by the appropriate factor 



,(b) 



/(I 



to arrive at 0.52 ± 0.24 ± 0.26 as the estimate of the contribution of events with 



one fake lepton to the total background. The uncertainties quoted are statistical and systematic, 
respectively, taking correlations into account. 

While the BTag-and-probe technique described above is based on different assumptions than 
the TL method of Section [5J| we note that in their implementation the two techniques are quite 
similar. Both techniques use events in a Rellso sideband combined with a scale factor deter- 
mined from an independent control sample to estimate the background in the signal region. 
The most notable differences are the requirement of the away-jet b-tag, which targets leptons 
from b decay in the BTag-and-probe method, the choice of variables used to parametrize ej/i 

in one case and e^j^jgg or e^^j^g in the other, and the size of the Rellso sideband used in the ex- 
trapolation. 

Combining the background estimates for events with one and two fake leptons and propagat- 
ing all statistical and systematic uncertainties between channels, including their correlations. 
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we arrive at a final estimate of background due to fake leptons of 0.70 ± 0.23 ± 0.21 events, with 
the first (second) uncertainty being statistical (systematic). 

5.3 Search using Hadronic Triggers and Xh 

Simulation studies show clearly that the largest source of background for the Zh channels is due 
to Th fake leptons. We estimate this background using the same "Tight-Loose" (TL) method as 



was used for fake leptons in Section 5.1 except that for the "Loose" selection we loosen the 
Th identification instead of the isolation. To be specific, part of the discrimination between 
hadronic T decays and generic QCD jets is based on five neural networks trained to identify 
different hadronic t decay modes. The neural network requirements are used for the tight, but 
not the loose selection. 



As in Section 5.1 in order to predict the number of events from fake r^, we measure the tight- 
to-loose ratio ejL in bins of rj and pr- On average, eji = 9.5 ± 0.5%, where the uncertainty 
is statistical only. This is measured using a single- T/, control sample with Hj > 300 GeV and 
£miss ^ 20 GeV. The Hj requirement results in hadronic activity similar to our signal region, 
while the E^^^^ requirement reduces contributions from W, Z plus jets, and tt, resulting in a 
control region that is dominated by QCD multijet production. 



Table 1: Validation of the TL method. The number of observed events is compared to the 
number of predicted events in simulation (first two columns) and in a background-dominated 
control region with relaxed selection criteria (last two columns). The simulation is normalized 
to 35 pb^^. The first and second uncertainties in the number of predicted events in data are 
statistical and systematic, respectively. 





Simulation 
Only SM 


Data 
Relaxed selection 


Channel 


Observed 


Predicted 


Observed 


Predicted 


TT 


0.08±0.03 


0.15±0.15 


14 


14.0±4.3±2.6 


eT 


0.35±0.12 


0.30±0.11 


1 


0.8±0.4±0.1 


f/T 


0.47±0.15 


0.49±0.20 


2 


2.9±0.6 ±0.4 



The expected number of background events is estimated by selecting er, }iT, and tt events 
where the t candidates pass the loose selection but fail the tight selection. We find 1, 2, and 2 
such events, respectively, in these three channels. For the eT and /yt channels, these events are 
weighted by the corresponding factors of cji/ (1 — eji^), while for the tt channel each event is 
weighted by the product of two such factors, one corresponding to each of the two t leptons in 
the event. 

We perform two types of validation of this background estimate. First, we compare observa- 
tion and prediction in the signal region for simulation. Second, we compare observation and 
prediction in data after relaxing the Hj selection from 350 GeV to 150 GeV, and removing the 
£miss requirement. In simulation the contribution of LMO represents less than 4% of events with 
two same-sign isolated leptons and Hj > 150 GeV. Table [ijpresents both of these validations of 
the background estimation technique. We find good agreement between the observation and 
prediction in all channels in data and in the simulation. 

The largest source of systematic uncertainties in the prediction of background events is due 
to lack of statistics of simulated events to validate the method (30%). In addition, we find 
uncertainties of 18%, 8% and 7% in T;,t;j, eXh, and }iTh, respectively, due to the correlation of Cjl 
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with Hj. We measure eji for Hj > 150 GeV. We determine the systematic uncertainty as the 
difference in the number of predicted background events from the reference measurement at 
Ht > 300 GeV with the measurement for Hj > 150 GeV. From simulation studies it is found 
that an additional 10% systematic uncertainty must be added in quadrature in the ez/, and /^(T;, 
channels to account for neglecting background contributions from fake electrons and muons. 
Taking all of this into account, we arrive at an estimate of 0.28 ± 0.14 ± 0.09 events with fake 
leptons, where the uncertainties are statistical and systematic, respectively. 

5.4 Comparison of Leptonic and Hr-Triggered Analyses with Electrons and 
Muons in the Final State 

As discussed above, we utilize two different trigger strategies to define selections that cover 
the maximum phase space possible in our search for new physics. In addition, the fact that 
these two selections have an overlap allows us to perform direct comparisons and cross-checks 
that we present in this section. 

We start by defining an overlap preselection requiring one electron or muon of pj >20 GeV 
and a second with pj > 10 GeV, Hj > 300 GeV, and no E^^^or isolation requirements. We note 



that this corresponds to the preselection sample from Section 5.2 with the pr of the leptons 



tightened to be consistent with Section 5.1 Comparing yields from this selection for the two 
trigger strategies on an event-by-event basis, we find that all of the H^-triggered events are 
also present in the lepton-triggered sample. From this comparison, we calculate efficiencies 
for the Hj and leptonic triggers to be (92 ± 4) % and ( 100^2 ) respectively. While statistics 
are limited, this confirms the trigger efficiencies measured in independent data samples, as 
presented in Section |4] 



Sections 5.1 and |5^ introduced two alternative methods for estimating the background due to 



fake leptons. Here we compare the two independent predictions in the region of overlap for the 
two searches. To compare the two methods in a common signal region, we require E^'^® > 30 



GeV and the lepton isolation described in Section 5.1 in addition to the preselection defined 
above. The TL method predicts 0.68 ± 0.39 based on a yield of 3 events that pass the loose 
selection. The second method introduced in Section [5^ results in an estimate of 0.27 ± 0.12 
events based on 11 events in the Rellso sideband. Both of these uncertainties are statistical 
only. We thus conclude that the two trigger strategies lead to consistent results within the 
kinematic region where they overlap, and the two methods of estimating backgrounds due to 
fake leptons give consistent results in that region. 

5.5 Electron Charge Mismeasurement 

A second potentially important source of background consists of opposite-sign dilepton events 
(e^e^ or e^f/^) where the sign of the charge of one of the electrons is mismeasured because of 
hard bremsstrahlung in the tracker volume. 

We measure the electron charge in three different ways. Two of the measurements are based 
on the reconstructed track from two separate tracking algorithms: the standard CMS track re- 
construction algorithm |29. 30| and the Gaussian Sum Filter algorithm [31 1, optimized for the 
measurement of electron tracks that radiate in the tracker material. The third measurement is 
based on the relative position of the calorimeter cluster and the projection to the calorimeter of 
a line segment built out of hits in the pixel detector. To reduce the effect of charge mismeasure- 
ments, we require agreement among the three measurements. 

After this requirement, the probability of mismeasuring the charge of an electron in simulation 
is at the level of a few per mille, even in the \ri\ > 1 region where the amount of material is 
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Figure 7: (Left) The probability to mismeasure the electron charge as a function of t] in the 
range 10—100 GeV, as obtained from simulation. (Right) Same-sign ee invariant mass distribu- 
tion in data compared with the Z — )■ ee expectation from simulation. 



largest, as can be seen in Fig. [7| (left). 

To demonstrate our understanding of this probability, we show in Fig. [7| (right) the invariant- 
mass spectrum for same-sign ee events and our simulation-based prediction from Z — > ee with 
one mismeasured charge. The Z sample shown uses the tight electron selection described in 
except with no jet or Hj requirement. Instead, we require £^'^® < 20 GeV and 



4.1 



Section 

transverse mass < 25 GeV to reduce backgrounds from W+ jets. The highest-pj lepton has 
been used in the calculation of the transverse mass. 

Measurement of the electron momentum is dominated by the energy measurements in the cal- 
orimeter, while the measurement of its charge is dominated by measurements in the tracker. An 
electron with mismeasured charge will thus still have a correctly measured momentum, lead- 
ing to the clear Z peak in the e^e"*^ invariant-mass displayed in Fig. [t] (right). Normalized to 35 
pb^^, the simulation predicts 7.4 ± 0.9 events in the Z mass region, with the quoted uncertainty 
due to statistics. In data, we observe 5 events in the same region. We predict the same-sign Z 
yield in simulation (6.37 ± 0.03) and data (4.9 ± 0.1) based on reweighting the opposite-sign 
Z — > e^e^ yield by a simulation-based parametrization of the probability for electron charge 
mismeasurement as a function of pr and //. 

For the leptonic trigger searches we estimate the number of background events due to charge 
misidentification by scaling the opposite-sign yields by the above probability function. We 
estimate the background due to electron charge misidentification as 0.012 ± 0.002 and 0.04 ± 
0.01 for the E^^'' > 80 GeV and Ht > 200 GeV regions, respectively 

For the Hr-triggered searches these backgrounds are further reduced since the opposite-sign 
yield is smaller given the tighter Hj requirement. The resulting background prediction is 
0.008 ± 0.005 (ee) and 0.004 ± 0.002 (eji) events. For the search with Tu in the final state, this 
background is negligible and ignored, as even in the opposite-sign er^ channel, background 
contributions dominate over those with a prompt T;,. 

We assign a 50% systematic uncertainty on the estimated backgrounds due to electron charge 
mismeasurement. This is motivated by the statistics available in the doubly charged Z — )■ e^e^ 
signal region. 
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6 Signal Acceptance and Efficiency Systematic Uncertainties 

Electron and muon identification efficiencies above pj ~ 20 GeV are known at the level of 3% 
per electron and 1.5% per muon, based on studies of large samples of Z — > ee and fifi events in 
data and simulation. The uncertainties increase as the efficiencies themselves decrease towards 
lower pt, reaching 6% (8%) per muon (electron) at 5 (10) GeV, based on studies of large samples 
of Z — )■ ee and in data and simulation. In addition, there is a potential mismodelling of 
the lepton isolation efficiency between data and simulation that grows with the amount of 
hadronic activity per event. To assess this, we compare the isolation efficiency as a function of 
track multiplicity in data and simulation for Z — )■ ee and mi, and extrapolate to new physics 
signals with large hadronic activity using simulation, as discussed in more detail in Section |8] 
Based on this, we assign an additional 5% systematic uncertainty per lepton. There is also a 1% 
(5%) uncertainty associated with the lepton {Hj) trigger efficiency. 

The efficiency of the hadronic t;, selection is studied in data via the process Z — > tt, where one 
T decays hadronically while the other decays into a muon [ISJ. The available statistics are an 
order of magnitude lower than the statistics available in Z — >• ee or fiji, at significantly lower 
purity. Accordingly, reconstruction versus pr can not be studied at the same level of detail 
in data as for the electron and muon reconstruction, and we depend to a greater extent on an 
accurate simulation than we do for electrons and muons. We assign an uncertainty of 30% [15] 
to the Tjj selection efficiency to account for the limited statistics available in data to validate the 
efficiency measured in simulation. 

An additional source of systematic uncertainty is associated with the current ~ 5% uncertainty 
on the hadronic energy scale Ill9| at CMS. This scale uncertainty limits our understanding of 
the efficiency of the Hj and E^'^® requirements. Clearly, final states where the typical Hj and 
E^^^^are large compared to the selection values used in the analysis are less affected than those 
with smaller Hj and £^^®^. We compute the systematic uncertainty due to this effect for the 
LMO benchmark point with the four signal selections using the method of Ref. 1(27} . We use the 
LMO model as it is typical of the possible SUSY final states to which these analyses are sensitive. 
We find that the uncertainty varies between 1% at Ht > 60 GeV and 7% at Ht > 350 GeV, the 
values of Hj used in the selections for the lepton-triggered baseline and t;, search respectively. 

Uncertainties in the acceptance due to the modelling of initial- and final-state radiation and 
knowledge of the parton density functions (PDF) are estimated to be 2%. For the latter, we use 
the CTEQ6.6 |32ll PDF and their uncertainties. 

Based on LMO as a signal model, we arrive at total uncertainties on signal efficiencies of 12%, 
15%, and 30% for the lepton triggered, Hj triggered low pr, and Hj triggered Th analyses, 
respectively. This includes a 4% luminosity systematic uncertainty p3|. In addition, to interpret 
these limits in terms of constraints on new physics models, one needs to take into account any 
model-dependent theoretical uncertainties. 

7 Summary of Results 

The results of our searches are summarized in Table |2] The background (BG) predictions are 
given by the rows labelled "predicted BG". In addition to the background estimate from data, 
we also present an estimate of the background based on simulation in the rows labeled as 
"MC". While QCD multijet production samples are used for testing background estimation 
methods in our control regions, they are too statistically limited to provide meaningful esti- 
mates of yields in the signal regions listed in Table |2j and are thus not included. All other SM 
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Table 2: Observed and estimated background yields for all analyses. The rows labeled "pre- 
dicted BG" refer to the sum of the data-driven estimates of the fake lepton contributions, and 
the residual contributions predicted by the simulation. The rows labeled "MC" refer to the 
background as predicted from the simulation alone. Rows labeled "observed" show the actual 
number of events seen in data. The last column (95% CL UL Yield) represents observed upper 
limits on event yields from new physics. 
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3.4 



simulation samples described in Section |4] are represented. Figure |8] summarizes the signal re- 
gion yields and background composition in all four search regions presented in Table |2] The 
lepton plus jets background where the second lepton candidate is a fake lepton from a jet clearly 
dominates all search regions. The low-p^-lepton analysis has a small, but non-negligible, back- 
ground contribution from events with two fake leptons. Estimates for backgrounds due to 
events with one or two fake leptons were obtained directly from data in appropriately chosen 



control regions, as described in detail in Sections 5.1 52} and 5.3 In the ee and e}i final states. 



small additional background constributions are present due to the electron charge mismeasure- 



ment, as discussed in Section 5.5 The remaining irreducible background from two prompt iso- 
lated same-sign leptons (WZ, ZZ, tfW, etc.) amounts to at most 10% of the total and is estimated 
based on theoretical cross section predictions and simulation. Uncertainties on the background 
prediction include statistical and systematic uncertainties added in quadrature. Contributions 
estimated with simulation are assigned a 50% systematic uncertainty. Data-driven estimates 
are assigned a systematic uncertainty between 30% and 50% across the various signal regions 
and channels. The ee,ef/, and channels have partially or fully correlated systematic uncer- 
tainties, as described in detail in Section |5] 

We see no evidence of an event yield in excess of the background prediction and set 95% CL 
upper limits (UL) on the number of observed events using a Bayesian method [34J with a flat 
prior on the signal strength and log-normal priors for efficiency and background uncertainties. 
These include uncertainties on the signal efficiency of 12% , 15%, and 30% for the lepton trig- 
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Figure 8: A visual summary of the observed number of data events, the expected number of 
background events, and the composition of the background for the four search regions. 



gered, Hj triggered low-pj, and Hj triggered T/, analyses, respectively, as discussed in more 
detail in Section |6j Based on the LMO benchmark model, the simulation predicts 7.3, 9.6, 9.1, 
and 2.0 events for the four signal regions, respectively. These LMO yields are based on the in- 
dividual NLO cross sections for all production processes that contribute to the expected signal 
yield. 

8 Interpretation of Results 

One of the challenges of signature-based searches is to convey information in a form that can be 
used to test a variety of specific physics models. In this section we present additional informa- 
tion that can be used to confront models of new physics in an approximate way by generator- 
level simulation studies that compare the expected number of events in 35 pb^^ with our upper 
limits shown in Table |2l 

The kinematic requirements described in Section |4] are the first key ingredients of such studies. 
The Hj variable can be approximated by defining it as the scalar sum of the pj of all final-state 
quarks (u, d, c, s, and b) and gluons with pj > 30 GeV produced in the hard-scattering process. 
The E^'^^can be defined as the magnitude of the vector sum of the transverse momentum over 
all non-interacting particles, e.g., neutrinos and LSP. The ratio of the mean detector responses 
for Hj and E^^^^as defined above, to their true values are 0.94 ± 0.05, and 0.95 ± 0.05, respec- 
tively, where the uncertainties are dominated by the jet energy scale uncertainty. The resolution 
on these two quantities differs for the different selections. In addition, the £^^^^ resolution de- 
pends on the total hadronic activity in the event. It ranges from about 7 to 25 GeV for events 
with Ht in the range of 60 to 350 GeV. The Hj resolution decreases from about 26% at 200 GeV 
to 19% for 300 GeV and to 18% for 350 GeV. The Hj resolution was measured in simulation 
using the LMO reference model, while the £™'®® resolution was measured in data. 

Figure |9] shows the efficiency versus pr using the LMO reference model for e,ji (left), and T;, 
(right). Efficiencies here include reconstruction, isolation, and selection. We fit the curves in 
Fig.|9]to the functional form: efficiency(pT) = e,nax + Ax (erf((pr - Pjcud/B) -1 ). We fix Pjcut 
to 10, 5, 15 GeV, and find {A, B, e^ax) of (0.40, 18, 0.66), (0.32, 18, 0.75), and (0.45, 31, 0.45) for e. 
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}i, and Tfi, respectively. 

Lepton isolation efficiencies depend on the hadronic activity in the event, and in some extreme 
cases like significantly boosted top quarks, on the event topology. The number of charged par- 
ticles at the generator level after fragmentation and hadronization is a good measure of the 
first of these two effects, and has been shown to agree reasonably well between data and sim- 
ulation II35I . We find that the isolation efficiency reduces roughly linearly by 10% for every 
15 charged particles with pj > 3 GeV within the detector acceptance. This was studied using 
Drell-Yan and LMO simulation, and compared with Z — )■ i^£^ data. The linearity was thus 
shown to be valid within a range of charged multiplicities from about 10 to 40. The LMO refer- 
ence model shown in Fig. |9]has an average charged multiplicity of ~ 25 for events that pass our 
selections. To arrive at the lepton efficiency for a new physics model with an average charged 
multiplicity of ~ 40, one would take the efficiency parametrization depicted in Fig.|9]and mul- 
tiply it by 0.9 to account for the 10% change in isolation efficiency due to the larger average 
charged multiplicity. The second effect, i.e., topologies as in boosted top quarks, is difficult to 
model at generator level, as the isolation efficiency may vary by an order of magnitude or more 
in extreme cases. Our results are thus not easily interpretable in new physics models with such 
characteristics without a detailed detector simulation. 
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Figure 9: Electron, muon (left) and Zh (right) selection efficiencies as a function of pj- The 
results of the fits described in the text are shown by the dotted lines. 



To show the level of precision obtainable based on such a simple efficiency model. Fig. 10 com- 
pares the exclusion line for CMSSM with tanjS = 3, Aq = GeV, and fi > 0, as obtained from 
simulation (solid blue) with the corresponding curve (dashed black) for this simple efficiency 
model. For the exclusion plane, the SUSY particle spectrum is calculated using SoftSUSY 1136 1 
along with sparticle decay using SDECAY II37I . The signal events are generated with PYTHIA 
6.4.22 [24J using CTEQ6m [38 1 PDF. The NLO cross sections are obtained with the cross section 
calculator Prospino L39J at each point in the {mo,mi/2) plane individually. 

When applying the simple efficiency model, we use the description above for the £^^^^ and 
Ht response and resolution, and the efficiency functions displayed in Fig. |9] No attempt was 
made here to correct for differences in hadronic event environment as those are small for the 
relevant {mo,mi/2) points. The width of the red shaded band around the blue line indicates 
uncertainties in the NLO cross section calculation. It includes variations of the PDF and simul- 
taneous variation by a factor of two of the renormalization and factorization scales. Both effects 
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are added in quadrature. Figure [To] shows that the theoretical uncertainties are larger than the 
imperfections in the simple efficiency model. The specific limit shown here corresponds to the 
leptonic trigger result with £^'^^ > 80 GeV for the purpose of illustration. The contour separates 
the bottom-left region where the expected event yield would be larger than the observed limit 
of 3.1 events (see Table|2| first row, last column) from the top-right region where such expected 
yield would be lower. 
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Figure 10: Exclusion contour in the mo — ^1/2 plane for CMSSM as described in the text. Com- 
paring the width of the red shaded band (theoretical uncertainty) around the blue curve with 
the difference between the solid blue and dashed black curves shows that the imperfections in 
the simple efficiency model described in the text are small compared to the theoretical uncer- 
tainties. 



We choose this CMSSM model because it provides a common reference point to compare with 
previously published Tevatron results ||40ll4l1l . The excluded regions from LEP are based on 
searches for sleptons and charginos Il42l - l46| . Other final states, especially the all-hadronic p2| , 
are better suited for this standard reference model, while the leptonic same-sign final state 
explored in this paper is more appropriate to constrain a wide variety of other new physics 
models HHZI- As discussed in Section|4j in a general supersymmetry context, one might expect 
gluino-gluino or gluino-squark production to lead to same-sign dilepton events via a decay 
chain involving a chargino. The salient, and very generic feature here is one lepton per gluino 
with either sign being equally likely. As a result, 50% of the dilepton events will be same-sign. 
Moreover, these cascade decays are typically characterized by two mass difference scales that 
separately determine typical Hj and lepton values. Different models of new physics may 
thus populate only one or the other of our different search regions. 
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9 Summary and Conclusions 

Using two different trigger strategies, we have searched for new physics with same-sign dilep- 
ton events in the ee, e}i, er^, fiTjj, and T/jT;, final states, and have seen no evidence for an ex- 
cess over the background prediction. The T;, leptons referred to here are reconstructed via their 
hadronic decays. The dominant background processes in all final states except T;,T;, involve 
events with one fake lepton. In the ZhTh final state, events with two fake Xh dominate. We have 
presented methods to derive background estimates from the data for all major background 
sources. We have set 95% CL upper limits on the number of signal events within \f]\ < 2.4 at 35 
pb^ in the range of 3.1 to 4.5 events, depending on signal region, and have presented details 
on signal efficiencies that can be used to confront a wide variety of models of new physics. Our 
analysis extends the region excluded by experiments at LEP and the Tevatron in the CMSSM 
model. 
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